en 



< 



International Journal of Theoretical Physics manuscript No. 

(will be inserted by the editor) 



Vector Models in VT Quantum Mechanics 



Katherine Jones-Smith • Rudolph 
Kalveks 



Received: 30 May 2012 / Accepted: 14 Jan 2013 



^-—^ Abstract We present two examples of non-Hermitian Hamiltonians which 

fSJ consist of an unperturbed part plus a perturbation that behaves like a vector, 

^ in the framework of VT quantum mechanics. The first example is a general- 

Q^ ization of the recent work by Bender and Kalveks, wherein the E2 algebra was 

examined; here we consider the E3 algebra representing a particle on a sphere, 
and identify the critical value of coupling constant which marks the transition 
fvq from real to imaginary eigenvalues. Next we analyze a model with S0(3) sym- 

metry, and in the process extend the application of the Wigner-Eckart theorem 
' _^ ' to a non-Hermitian setting. 

t^ Keywords Non-Hermitian quantum mechanics • PT quantum mechanics • 

+^ Wigner-Eckhart theorem 

a PACS ll.30.Er • 03.65.-w • 02.30.Fn 

1 Introduction 

i^ There are many situations in quantum mechanics wherein the Hamiltonian 

i^ under consideration can be written as 

\o 
in 



H = Ho + H, (1) 
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where Hq is the unperturbed part and commutes with the generators Ti of 
symmetry group G: 

[Ho,T,] = (2) 

and Hi can be treated hke a perturbation and behaves hke a vector under G. 
We wish to generahze this situation in the context of VT quantum mechan- 
ics [HE], where the assumption that operators such as the Hamiltonian are 
Hermitian is relaxed, and replaced by other requirements, notably that the 
Hamiltonian commutes with the parity {V) and time-reversal (7") operators. 

Interest in non-Hermitian quantum mechanics continues to grow [5], and 
recently a number of experiments have observed the so-called VT phase tran- 
sition, where the eigenvalues of a VT Hamiltonian make a transition from 
being complex to real once a critical value of a coupling constant is reached 
[1], [S], [B]. Thus it is relevant to seek new T^T-counterparts to conventional 
Hamiltonians. 

In this work we present two simple cases that can be described as non- 
Hermitian vector perturbation models where the Hamiltonian can be written 
as in eq (fl]); first we consider a particle confined to the surface of a sphere, 
where the Hamiltonian acts within an infinite dimensional Hilbert space, and 
next we consider a generic vector perturbation within a finite dimensional 
Hilbert space and determine the spectrum of eigenvalues using the Wigner- 
Eckart theorem. We find that for a range of parameters each of these models 
has a pure real spectrum. At critical values of the coupling the model undergoes 
VT transitions wherein the eigenvalues become complex. 



2 E3 Algebra: particle on a sphere 

We begin by generalizing the analysis presented in [7]. They considered the 
E2 algebra which consists of elements J, u, v such that 

[J, u] = —iv [J, v] = iu [u, v] = 0. (3) 

The Hamiltonian 

h = J'^ + igu (4) 

where J — —id/ 89, u — sin0, v — cos 6 and g is a constant, represents a 
2-dimensional quantum particle restricted to radius r — 1. 

A generalization of this is the E3 algebra and restricting the particle to the 
surface of a sphere (r — 1). This is described by the Hamiltonian 

h = L^ + iguz (5) 

where L obeys 

[Li,L.j] = ieijkLk (6) 
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M is a vector operator whose components are given by 

(7) 
(8) 
(9) 

and 5 is a constant. The remaining commutators are straightforward to calcu- 
late; 

[Li,u.j\=ieijkUk K,Uj]=0. (10) 

Following Bender and Kalveks we consider the case of even time reversal: 
for a wave function ip {9,(f>) the time reversal operator T is manifested as 
complex conjugation: 

T^{e,<j,)^r{o,<l>) (11) 

hence T^ = 1. It is easy to verify the action of T on the elements of the 
algebra: TLfT = —Li and Tuff — Ui. The parity operator V takes ^ into the 
antipodal point: 

Pi:{e,<t))^ip{-K-e,(i) + TT) (12) 

so V^ = 1; elements transform under parity as VLiP = Li and VuiV = —Ui. 
Note that the Hamiltonian ft, in eq ([5]) commutes with the combined operation 
VT but not with V or T individually. Now let us determine the eigenvalue 
spectrum of this Hamiltonian. We wish to solve 

hij{e,cp) = x^{0,<f>) (13) 

and we try the general solution: 

^ ((?,</>) = / (0)6™^. (14) 

For convenience we define ry — cos9; this simplifies the eigenvalue equation for 
/: 

- (1 - ^') a^ + 2^7^ + Y^,f + WVf - A/ (15) 

where ttt. is a fixed integer. If we let 

(16) 





ho = 


-(1 


-"''0-'|- 


m2 




l-ry2 


then the Hamiltt 


Dnian we 


wish to solve is 












hnf + igr]f = 


A/. 





(17) 

We impose the boundary condition that the solution must be regular at 77 = 
±1. 

Let us choose basis elements 

\l) ^ NiP^„,\{7^) (18) 
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Fig. 1 Real and imaginary components of eigenvalues E for the Hamiltonian given by 
Eq |17| The first six eigenvalues for m = (blue) and m = 1 (green) are shown. Intercepts 
on the E axis are given by l!(£ + 1) for i! = to 6. For the case of m = 0, we find that the 
spectrum is entirely real for < g < 1.899 at which point there is a transition to one pair of 
complex conjugate eigenvalues in the spectrum. At <; = 11.45 there is a second transition, 
to two pairs of complex conjugate eigenvalues. Similarly for the case of m = 1, we find 
one complex conjugate eigenvalue pair at g = 5.41, and two pairs at g = 19.04. In these 
computations the Hamiltonian is truncated to a 100x100 matrix; we have verified that the 
relevant part of the spectrum is insensitive to the truncation. 



where I 



\m\, \m\ + 1, 



Pi \m\ are the associated Legendre polynomials, 



with conventional normalization factor 



N,^ 



(21 + 1) {i^\m\y- 

2 y{l + \m\)l- 



(19) 



hQPL\m\iv)^lil + '^)PL\m\{v) 



(20) 



The Plum's satisfy 

SO the matrix of ho in this basis is diagonal. The matrix elements of the poten- 
tial term, igi], can easily be determined from the normalization and recursion 
relations of the P/^^'s. By diagonalizing the truncated Hamiltonian matrix we 
can numerically obtain the eigenvalues of eq (17 1; see Fig. (111). 



3 T'T" Vector Model in Finite-Dimensional Hilbert Space 



E3 may also be regarded as a realization of the VT vector model with symme- 
try group 5*0(3) and for which the Hilbert space is infinite dimensional. Now 
we wish to turn out attention to realizations of the VT vector model with 
finite dimensional Hilbert spaces. Let us write a simple, generic Hamiltonian 
H = Ho + Hi where 



Hf) — L„ 



Lt 



Hj = V, 



(21) 
(22) 



and Vz is the z component of a vector operator. 

Our task is to obtain a matrix representation of the total Hamiltonian, solve 
for its eigenvalues and determine what value of the non-Hermitian perturbation 
cause the eigenvalues to become complex. 
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Naturally we choose to work with the angular niomentum eigenstates 
1^, m); the action of Hq on these states is well known, and we can utilize 
the Wigner-Eckart theorem to determine the action of Hi = V^ ■ 

Note that the dimensionality of the relevant vector space depends on the 
angular momenta of the multiplcts but clearly it is finite. Suppose we consider 
the two multiplcts \£,m) and \(, + 1,to); to takes on values from —t to +(, in 
the first multiplet and from — £ — 1 to £ + 1 in the second multiplet, so there 
are (2£ + 1) + (2£ + 3) = 4^ + 4 of these states. 

The action of Hq on these states is simply 

L^\i,m) ^i{l + l)\i,m) (23) 

L2|^+1,to) ==(^ + l)(£ + 2)|^+l,m). (24) 

So all that remains is to determine how Vz acts on these states; here we employ 
the Wigner-Eckart theorem, which we have extended to the non-Hermitian 
case as detailed in Appendix A. We find (f',TO-'|y2|€, to) = unless to = m! . 
Thus we need only to determine 

(^,to|v;|^,to), 

(^ + 1,to|14|^ + 1,to), 

{t,m\Vz\t+ 1,to), and 

(^ + 1,to|K|^,to) 



in order to completely specify V^ in this space. The first two in this list can 
be expressed in terms of the reduced matrix element a defined in Appendix 
A; in general we find 



(€, to-IFjI^^jto) — mai 
(^+1,to|V;|^+1,to) =ma2; 



however we also wish to enforce WzV — —Vz and TVzT — —Vz-, which 
restricts oi\ = ot^ = 0. (Determination of V and 1' within this space follows 
straightforwardly from their action on the spherical harmonics PYimiO,4>) — 
{-lYY,^i9,<l>) and TYimie^) = Y;^i0,<l>) = (-l)™r,,_™(0, ,/)).) 

For the other two types of matrix elements, {£,m\Vz\i + 1,to) and(^ + 
1, to-JFjI^, TO.), we find these are proportional to other reduced matrix elements 
/3 and 7; 

{e+l,m\Vz\e,m)^fiml3, 



where 

(^+1)2 -to2 T^/" 

Jim — 



(2£+l)(2£ + 2) 



(25) 
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Note that fim is even in m. When we enforce PVzP ^ —Vz and TVzT = — T4, 
we find this requires /? and 7 to be pure imaginary, so we define l3 = ib, j = ic 
for some real numbers 6, c. Note that in determining these matrix elements we 
do not assume V is Hermitian; we rely only on the commutators of V with the 
angular momentum operators. (See Appendix A for details.) 

Now let us write down the matrix corresponding to the Hamiltonian H = 
Ho + Hj. Consider the two-dimensional subspace spanned by the states \£,m) 
and \£ + l,m) for a fixed value of m that lies in the range —£, ...,£. Within 
this subspace 

[£{£ + 1) 

(^ {£+!){£ + 2) 



Hn^ 



(26) 



and 



Vz 





-ibfim 







(27) 



In addition consider the two dimensional subspace spanned by the states \£-\- 
1, £ + 1) and 1^+1,-^—1). These states are not coupled by the perturbation 
Vz to any other state and hence Vz = within this subspace. On the other 
hand the unperturbed Hamiltonian in this subspace is given by 



Hn = 



(^+l)(^ + 2) 

{£ + l){£ + 2) 



(28) 



It is convenient to define 



hi+i — 



{£ + l){£ + 2) 






(^+l)(^ + 2) 



(29) 



and 



where m = 
matrix 



f^m. — 



£{£ + 1) -icfir^ 



(30) 



-ibfi,^ {£ +!){£ + 2)^ 
—£, ...,£. The Hamiltonian can now be written as a block-diagonal 



( hi+i 



V 



hi-i 



\ 



'J 



(31) 



The individual 2x2 matrices that constitute the Hamiltonian are sim- 
ple enough that we can obtain analytic expressions for the eigenvalues. The 
eigenvalues of ht+i are two-fold degenerate and arc simply (£+ V){£ + 2). The 
eigenvalues of hm are 



Xf^ = {£+iY±\{£+i?-hcf, 



(32) 
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Fig. 2 Real and imaginary parts of the eigenvalues Ai^^ assuming b = c. The blue line at 
A = 6 corresponds to the 2x2 matrix denoted hg^i in the text, with eigenvalues (£+l)(£+2). 
The other eigenvalues are m-dependent and correspond to the 2x2 matrices denoted hm 
in the text, with eigenvalues given by eql|32|l. As noted A^ „j = Af _„, so there are only two 
distinct m-dependent curves for i = 1. In each figure, \m\ = 1 is plotted in red and m = 
is plotted in green. Note that the transition to complex eigenvalues occurs at 6 = 4 for \^ ^ 

and b - 



12 Si 3.47 for Af g. 



Note that X^^m = Xi-m so for all m 7^ the eigenvalues of hm are also two-fold 
degenerate. Clearly, the eigenvalues are real provided 



be < 



■' Irn 



(33) 



We can make the following observations about the behavior of the eigenvalues. 
Once VT symmetry is broken, A^ and A~ form a complex conjugate pair. Since 
fim has its maximum value for m = 0, X^ becomes complex for m = first. 
Similarly, f^m is minimal for \m\ = £, so A^ so these are the last eigenvalues 
to go complex. For example we consider the case £ = 1 and choose b = c for 
simplicity. We plot the eigenvalues in Figure [2] 

It is worth noting that in the Hermitian case b ~ — c. Hence the condition 



in cq (33 1 that ensures the eigenvalues are real is always met. 



4 Conclusion 



We conclude by noting two natural generalizations of our results that deserve 
further investigation. First the model of a particle on an ordinary 2-sphere 
considered in section II may be generalized to a particle on a sphere in n 
dimensions. The VT transition for this model may be amenable to analytic 
study in the large n limit and may shed some light on VT symmetric non- 
linear sigma models of which it would represent a -)- 1 dimensional case [8] . 
Second the vector model constructed in section III may be easily generalized 
from the symmetry group S0(3) to any Lie group and therefore represents 
only one member of a large class of such models. 
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Appendix A: Wigner-Eckart Theorem 

Suppose we have an angular momentum operator L and a vector operator 
V satisfying the commutation relations 



[Li,Vj] = ieijkVk- 



(34) 



Let \i, m) denote an angular momentum multiplet of total angular momentum 
f. and z-component m. Then according to the Wigner-Eckart theorem the ma- 
trix elements of Vz and V± = VxUVy between multiplet states are determined 



by the commutation relations eq ( 34 1 . In the usual Wigner-Eckart theorem the 



Cartesian components of the operator V are assumed to be hermitian. Here 
we present a non-Hermitian generalization of the theorem. 
Following the usual arguments wc find the selection rules 



{£' ,m'\Vz\£,m) = unless m! ^ m 
(£', m'\Vj^\t, m) = unless m' = to + 1 
(^', m'\V- \i, m) — unless m' — m— 1. 



(35) 
(36) 
(37) 



Furthermore the matrix elements vanish unless I' = i—1 or £' = i or i' = i+1. 
Consider the case £' = £. Generalization of the usual arguments shows that 

{£, m + l\V+\i, to) = A{e - m)^/^{e + m + 1)^/^ m = -£,... J - 1 

{e,m\Vz\i,m) = Am m^ -£,...,£ 

{i,m- l\V^\i,m) = A(£-to + 1)1/2(^ + to)1/2 m = -^+1,...,^ 

(38) 

where the proportionality constant A is a complex number called the "reduced 
matrix element" . Note that for V hermitian, A would have to be real, but there 
is no such restriction in the non-hermitian case. 
Similarly in the case ^' = ^ -f 1 we find 



{£ + l,m + l\V+\£,m) ^ B 

{£+l,m\Vz\£,m) ^ ^B 
{£+l,m-l\V-\£,m) = -B 



{l + m + 2){l + m + iy 


(2£ + 2)(2^-f 1) 


' (e - m + l){£ + m + 1) 


{2£ + 2){2i+l) 


'{e-m + l){£-m + 2) 



1/2 



(2^ + 2)(2£+l) 



1/2 



1/2 



(39) 



where m 



-t 



and B is another complex reduced matrix element. 
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Finally in the case that I' 



1 we find 



{(.~^l,m + l\V+\(.,m) =-C 

(£-l,m|T4|£,m) = -C 

{t-l,m-\\V-\l,m) = C 



{I - m - l){i - m)' 
(2l){2t-l) 



1/2 



{2l){2t-l) 

{i + m){i + m- 1) 
{2l){2l-l) 



-,1/2 



(40) 

where C is a complex reduced matrix element and m ~ —£, ...,£ — 2 in the 
first line of eq (401, m — ~£ + 1, ...,£— 1 in the second line of eq (40), and 
m = —£ + 2, . . . ,£ in the last line of eq ( 40 1 . 



In the hermitian case the reduced matrix elements satisfy B = C* but in 
the non-hermitian case there is no such restriction on the complex elements B 
and C. 
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